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Abstract 



We present the results of theoretical investigation on the electronic struc- 
ture, bonding nature and ground state properties of Tl^Al and TI12AIH4 using 
generalized-gradient-corrected first-principles full-potential density-functional 
calculations. TI12AIH4 has been reported to violate the "2 Arule" of H-H sep- 
aration in hydrides. From our total energy as well as force-minimization 
calculations, we found a shortest H-H separation of 1.95 A in accordance with 
recent high resolution powder neutron diffraction experiments. When the 
TI12AI matrix is hydrogenated, the volume expansion is highly anisotropic, 
which is quite opposite to other hydrides having the same crystal structure. 
The bonding nature of these materials are analyzed from the density of states, 
crystal-orbital Hamiltonian population and valence-charge-density analyses. 
Our calculation predicts different nature of bonding for the H atoms along a 
and c. The strongest bonding in TI12AIH4 is between Th and H along c which 
form dumb-bell shaped H-Th-H subunits. Due to this strong covalent interac- 
tion there is very small amount of electrons present between H atoms along c 
which makes repulsive interaction between the H atoms smaller and this is the 
precise reason why the 2 A rule is violated. The large difference in the inter- 
atomic distances between the interstitial region where one can accommodate 
H in the ac and ab planes along with the strong covalent interaction between 
Th and H are the main reasons for highly anisotropic volume expansion on 
hydrogenation of TI12AI. 

PACS numbers: 71., 81.05.Je, 71.15.Nc, 71.20.-b 
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I. INTRODUCTION 



Hydrides of intermetallics have been extensively studied because of their applications in 
re-chargeable batteries. Unfortunately, most metals that absorb large amounts of hydrogen 
are either heavy or expensive.0 Consequently, there is a constant search for hydrides that 
may be suitable for practical applications. First of all, it is very important to understand 
how crystal structural evolution takes place in the course of hydrogenation. Numerous 
studies have been done to explain observed stabilities, stoichiometries, and preferred H sites 
in hydrides of metallic and intermetallic compounds. Structural studies of hydrides have 
provided empirical rules! that can be used to predict the stability of the H sublattice in a 
given metal configuration. A survey of stable hydrides show that the H— H distance does 
not go below 2.1 A (the 2 A rule) with a minimum radius of 0.4 A for the inter-site to be 
used for the accommodation of H. These rules have been used to predict new hydrides whose 
existence is verified experimentally.liH! 

The review of Yvon and Fischer! states that Th 2 AlH 4 l and K 2 ReHglil are violate the 
2 A rule, the shortest H— H separation being 1.79, and 1.87, respectively. K 2 ReHg is classi- 
fied among complex transition metal hydrides, which comprise highly covalent solids with 
nonmetallic properties. TI12AIH4, on the other hand, has metallic character. 

ThoAlS together with Zr 2 Fe, Zr 2 Co and Zr 2 Ni crystallize in the CuAl 2 -type structure, 
whereas their hydrides form rather different crystal structures. Zr 2 Fe and Zr 2 Co form the 
isostructural deuterides Zr 2 MD 5 (M = Fe,Co)@ with a change in symmetry from 1 "4/ 'mem 
to PA/ncc on deuteration. Th 2 AlH_4 and Zr 2 NiH-4.740 are formed without any change in 
the symmetry from their parent structures. Th 2 AlH4 belongs to the exclusive class which 
does not obey the 2 A rule. The lattice expansion along a and c has proved to be highly 
anisotropic on hydrogenation of Th 2 Al. In order to shed light on this effect we need theo- 
retical understanding about bonding nature in this compound. Further, the understanding 
of the lattice expansion and distortion during hydrogenation will be important for the eval- 
uation of stability of the hydride. So, we have made detailed study of Th 2 Al and Th 2 AlH4 
by first-principle calculations. 

Two different powder neutron diffraction (PND) of Th 2 AlH 4 give different H-H 
separations, viz. the older! value is 1.79 A and the more recent! value is 1.97 A. So one 
aim of this study has been to solve this discrepancy. In principle, the stability of hydrides 
can be evaluated directly from a theoretical study of the total energy. However, owing to 
the complexity of the structure of transition metal hydrides, no reliable theoretical heat of 
formation has hitherto been reported.0 Nakamura et a/.0 were the first to calculate heat 
of formation. However, these authors obtained positive and unrealistically large heat of for- 
mation even for stable La-Ni based hydrides except for (La 2 Ni 10 H 14 ).liaEJ This unfavorable 
result clearly indicates that local relaxation of the metal atoms surrounding the hydrogens 
must be included in the calculations in order to predict the structural stability parameters. 
Hence our calculations take into account local relaxation by optimizing the atom positions 
globally. 

We present the electronic structure of Th 2 Al and Th 2 AlH 4 , obtained by the full-potential 
linearized-augmented plane wave (FPLAPW) method. A central feature of the paper is the 
evaluation of the electronic structure and bonding characteristics on introduction of H into 
the Th 2 Al matrix. In addition to regular band-structure data, we also provide crystal 
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orbital Hamiltonian population (COHP)0E1 results to illustrate the chemical bonding in 
more detail. 

This paper is organized as follows. Details about the involved structure and computa- 
tional method are described in Sec. ||. Sec. |TTT| gives the results of the calculations and 
comparisons with the experimental findings. Conclusions are briefly summarized in Sec. [TV]. 



II. STRUCTURAL DETAILS 

Th 2 Al and Th 2 AlIFi crystallize in space group 14 /mem with the lattice parameters a 
= 7.618, c = 5.862 A0 for Th 2 Al and a = 7.626,c = 6.515 A for T1i 2 A1hJ The crystal 
structure of Th 2 AlH 4 is illustrated in Fig. 0. The crystal structure of Th 2 Al contains four 
crystallographically different interstitial sites, which are the suitable sites for hydrogen ac- 
commodation, 16/ and 46 each coordinated to four Th, 32m coordinate to three Th and 
one Al 16k coordinate to two Th and two Al. Each 16/ based intersite tetrahedra share a 
common face with another 16/-based tetrahedron, whereas the 46-based tetrahedron shares 
each of its four faces with 16/-based tetrahedra. Some of the tetrahedral intersites are closely 
separated owing to the face sharing of the coordination polyhedra. According to the exper- 
imental findings,!'! the 16/ sites are fully occupied in Th 2 AlD 4 , and also the structure is 
completely ordered. 



A. Computational details 

In our calculations we use the full-potential linearized augmented plane wave (FP- 
LAPW) method in a scalar relativistic version without spin-orbit coupling as embodied in 
the WIEN97 code.0 In brief, this is an implementation of density-functional theory (DFT) 
with different possible approximations for the exchange and correlation potential, including 
the generalized-gradient approximation (GGA). The Kohn-Sham equations are solved using 
a basis of linearized augmented plane wavesJll For the exchange and correlation potential, we 
used the Perdew and Wangil implementation of GGA. For the potential and charge density 
representations, inside the muffin-tin spheres the wave function are expanded in spherical 
harmonics with l max = 10, and non spherical components of the density and potential are 
included up to l max = 6. In the interstitial region they are represented by Fourier series and 
thus they are completely general so that such a scheme is termed full-potential calculation. 
The present calculations we used muffin-tin radii of 2.5, 2.1 and 0.9 Bohr for Th, Al and H 
respectively. 

The basis set includes 7s, 7p, 6d and 5/ valence and 6s and 6p semi-core states for Th, 
3s, 3p valence and 2p semi-core states for Al and Is states for H. These basis functions 
were supplemented with local orbitals0 for additional flexibility to the representation of the 
semi-core states and for generalization of the linearization errors. We have included the local 
orbitals for Th 6s, 6p and Al 2p semicore states. In all our calculations we have used tetra- 
hedron method on a grid of 102 k points in the irreducible part of the hexagonal Brillouin 
zone (IBZ),@ which corresponds to 1000 k points in the whole Brillouin zone. The calcula- 
tions are done at several cell volumes (around the equilibrium volume) for both Th 2 Al and 
Th 2 AlH 4 and corresponding total energies are evaluated self-consistent ly by iteration to an 
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accuracy of 10 -6 Ry./cell. Similar densities of k points were used for the force minimization 
and c/a optimization calculations. 

In order to measure the bond strengths we have computed the COHP0 which is adopted 
in the TBLMTO-47 packageSEI COHP is the density of states weighted by the corre- 
sponding Hamiltonian matrix elements, which if negative indicates a bonding character and 
positive indicates an anti-bonding character. The simplest way to investigate the bonding 
between two interacting atoms in the solid would be to look at the complete COHP between 
them, taking all valence orbitals into account. However, it may sometimes be useful to focus 
on pair contributions of some specific orbitals. 



III. RESULT AND DISCUSSION 



The H— H separation is one of the most important factors in identifying the potential 
candidate for hydrogen storage, because if the H-H separation is small one can accommodate 
more H within a small region. From this point of view,Tli2AlH4 may be considered as a 
potential candidate for storing H. To the best of our knowledge no theoretical or experimental 
attempts have been made to study cohesive properties like heat of formation (AH), cohesive 
energy (E coh ), bulk modulus (B ) and its pressure derivative (B' Q ) for this compound. Hence, 
this is the first theoretical attempt to study the ground state properties and bonding in this 
compound. 



A. Structural optimization from total energy studies 



In order to analyze the effect of hydrogenation on the crystal structure of Th 2 Al and 
to verify the discrepancy between the experimentally observed H-H separation, we have 
optimized the structural parameters for TI12AI and TI12AIH4. For this purpose, first we have 
relaxed the atomic positions globally using the force-minimization technique, by keeping 
experimental c/a and cell volume (Vo) fixed to experimental values. Then the theoretical 
equilibrium volume is determined by fixing optimized atomic positions and experimental cj a, 
and varying the cell volume by ±10 % of Vq. Finally the optimized c/a ratio is obtained by 
a ±2 % variation in c/a ratio (in steps of 0.005), while keeping the theoretical equilibrium 
volume fixed. It is important to note that experimentally observed lattice parameters are 
almost same, while the atomic position of H alone differs between the two experimental 
results (according to Bergsma et aM H coordinates are 0.368, 0.868, 0.137 and Sorby et a/.i 
give 0.3707 , 0.8707, 0.1512). The total energy vs. cell volume and c/a ratio curves for 
TI12 Al and TI12 AIH4 are shown in Figs. ||] and || respectively. From these illustrations it is 
clear that the equilibrium structural parameters obtained from our theoretical calculations 
are in very good agreement with those obtained recently by PND measurement s.@ 

The optimized atomic positions along with the corresponding experimental values are 
given in Table|I[ Table[H] gives calculated lattice parameters and interatomic distances, along 
with corresponding experimental values for both TI12AI and TI12AIH4. The theoretically 
estimated equilibrium volume is underestimated by 0.27% for TI12AI and 1.8 % for TI12AIH4. 
The underestimation of bond length in the present study is partly due to the limitation of 
local density approximation used in the calculations and also neglect of the zero-point motion 
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and thermal expansions. The difference between the experimental values may be due to the 
poor resolution of the earlier(1961) PND datJ. 

B. Cohesive properties 

The method of calculation for cohesive properties for intermetallic compounds are well 
described in Refs. Effiffil The cohesive energy is a measure of the force that binds atoms 
together in the solid state. The cohesive energy of a system is defined as the sum of the total 
energy of constituent atoms at infinite separation minus the total energy of the particular 
system. This is a fundamental property which has long been the subject of theoretical 
approaches. The chemical bonding in intermetallic compounds is a mixture of covalent, 
ionic and metallic bonding and therefore the cohesive energy can not be determined reliably 
from simple models. Thus, first principle calculations based on DFT have become a useful 
tool to determine the cohesive energy of solids. For the study of phase equilibrium the 
cohesive energy is more descriptive than the total energy, since the latter includes a large 
contribution from electronic states that do not play a role in bonding. From our cohesive 
energy calculations we get E co h = 0.15, 0.185 eV/atom for Th 2 Al and Th 2 AlH 4 respectively, 
indicating that hydrogenation enhances the bond strength in Th 2 Al. 

The formation energy (AH) is introduced in order to facilitate a comparison of system 
stability. AH is defined as the total energy difference between the compound and weighted 
sum of the corresponding total energy of the constituents. For the AH calculations, we 
used the total energy value of 2.320 Ry for the hydrogen molecule which was calculated 
with the von Barth-Hedin exchange- correlation potential.^ AH provides information about 
the stability of Th 2 Al towards hydrogenation. The calculated AH values for La-Ni based 
hydrides,0 were almost double the experimental^^ AH. As the LMTO-ASA method was 
used in this study, this discrepancy is expected because the internal relaxation of the atoms 
was not taken into account and the interstitial potential is not well represented in LMTO- 
ASA method. Therefore, AH calculated by using the full-potential method should be more 
reliable. Our calculated values for E co h and AH are given in Table fH\. Since (AH) is more 
negative and E co h is higher for Th 2 AlH 4 than for Th 2 Al, we can conclude that Th 2 AlH4 is 
more stable than Th 2 Al. However, no experimental AH values for Th 2 Al and Th 2 AlH4, are 
available, but we note that our calculated AH is close to the experimentally observed values 
of other Th-based hydrides, like ThNi 5 H 4 , TI1C05H4 and ThFeshLi all having the AH value 
of -36.63 kJ/(mol-H).i 

The bulk modulus of Th 2 Al and Th 2 AlH 4 was obtained by self-consistent total energy 
calculations for 8 different volumes within the range of V/Vq from 0.75 to 1.10 (see Fig. 
From the derivative of total energy with respect to the volume, the calculated bulk mod- 
ulus for Th 2 Al is 93.42 GPa and for Th 2 AlH 4 is 111.36 GPa. The corresponding pressure 
derivative of the bulk modulus (B' Q ) are 3.43 and 3.51, respectively. The enhancement in B 
value in the hydrogenated phase indicates that hydrogen plays an important role in bonding 
behavior of Th 2 AlH 4 . In particular, the hydrogenation enhances the bond strength, and 
hence the change in volume with hydrostatic pressure decreases with hydrogenation. This 
conclusion is consistent with the observation made from our calculated heat of formation 
and cohesive energy for Th 2 Al and Th 2 AlH 4 . 
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C. Anisotropic behavior 

For compounds which maintains the basic structural frame work, the occupancy of hy- 
drogen in interstitial sites is determined by its chemical environment (different chemical 
affinity for the elements in the coordination sphere also results in different occupancy). Al- 
though the H atom is small and becomes even smaller by chemical bonding to the host, it 
may deform and stress the host metal considerably depending upon the chemical environ- 
ment. Lattice expansion usually of the order of 5 to 30%, often anisotropic, results from 
hydride formation. The maximum volume expansion observed for CeRu2 to CeR^Ds (37%) 
is due to a hydrogen induced electron transition as shown by XPS measurements £21 A lat- 
tice contraction upon hydrogenation has so far only been observed for ThNi2 to ThNi2-D2 
(2.2%). For most hydrides formed from intermetallic compounds the crystal structure usu- 
ally changes with a loss of symmetry.0 In general the symmetry decreases as a function of 
hydrogen content and increases as a function of temperature. However, on hydrogenation 
of Th 2 Al the symmetry remains unchanged. 

The volume expansion during hydrogenation of Th 2 Al is 12.47% (AV/H atom is 
10.32A 3 ). This volume expansion is strongly anisotropic and proceeds predominantly per- 
pendicular to the basal plane of the tetragonal unit cell; Aa/a = 0.026%, Ac/c = 12.41%. 
This indicates a relatively flexible atomic arrangement in the [001] direction. In spite 
of the isostructurality between Th 2 Al, Zr 2 Fe (hydrated: Zr 2 FeH 5 )i and Zr 2 Co (hydrated: 
Zr 2 CoH 4 . 82 )t-3 the latter two exhibit a quite opposite anisotropic behavior in that the the 
unit cell expands exclusively along the basal plane. The c/a ratio plays an important role for 
the structural properties of intermetallic compounds including metal hydrides. For example, 
in the case of Zr 2 Fei, Zr 2 Coll, Zr 2 Ni0>fl and Th 2 Al c/a is 0.878, 0.867. 0.812 and 0.7695 
respectively, and for the corresponding hydrides Zr 2 FeH 5 ,0 Zr 2 CoH 4 §2," Zr 2 NiH 4 74 and 
Th 2 AlH 4 c/a is 0.810, 0.815, 0.833 and 0.8543 respectively (see Fig.gj). The increase in c/a 
for Zr2CoH 4 .g2 and Z^FeHs compared with the corresponding unhydrated parents is smaller 
than that for other compounds. On hydrogenation, the increase in c/a ratio for TI12AI is 
considerably larger than for Z^Ni, which may be the reason why the former retains the 
symmetry on hydrogenation. Our calculations describe well the anisotropic changes in the 
crystal structure on hydrogenation of Th 2 Al (see Table||). The c/a ratio increases almost 
linearly (Fig. ^) on going from Z^Fe to TI12AI whereas the corresponding hydrides show the 
opposite behavior. Hence, it appears that the systematic variation in c/a plays a major 
role in deciding the crystal structure for the CuAl2-type hydrides. When c/a < 0.825 the 
symmetry is changed from 1 '4/ 'mem to PA/ncc on hydrogenation, whereas when c/a > 0.825 
the crystal symmetry is apparently not affected. 

D. Electronic structure 

In order to understand the changes in the electronic bands on hydrogenation of Th 2 Al 
we show the energy band structure for Th 2 Al and Th 2 AlH 4 in Fig.|5]a and b respectively. 
The illustration clearly indicates that inclusion of H in the TI12AI matrix has a noticeable 
impact on the band structure, mainly in the valence band. The two lowest lying broad 
bands in Fig. |5|a originate from Al 3s electrons. As the unit cell contains two formula units, 
eight electrons are additionally introduced when Th 2 AlH 4 is formed from Th 2 Al. These 
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electrons form four additional bands (Fig.[5|b), a large deformation of the band structure 
is introduced by the of hydrogen in the Th 2 Al lattice. These bands become localized, the 
lowest lying energy band is moved from —7.34 to— 9.20 eV, and, the character of the latter 
band is changed from Al 3s to H Is character. The Al 3s bands are located in a wide 
energy range from —2.8 to — 7.34eV in Th 2 Al and are in a narrow energy range from —2.5 
to — 4.2eV in the hydride. The drastic change in the Al bands on hydrogenation of Th 2 Al is 
due to the electron transfer from H to Al and this is discussed further in Sec. |111 F| . The H s 
bands are found in the energy range from —2 eV to the bottom of the valence band. Their 
contribution at Ep is negligibly small indicating the formation of localized bands. The bands 
at Ep is dominated by the Al 3p and Th 6d electrons in both Th 2 Al and Th 2 AlH 4 . Owing 
to the creation of the pseudogap feature near Ep, the contribution of the Al 3p electrons to 
the bands at Ep level are significantly reduced by the hydrogenation of Th 2 Al. 



E. Nature of Chemical Bonding 

1. Density of state 

In order to obtain a deeper insight into the changes in chemical bonding behavior on 
hydrogenation of Th 2 Al we give the angular-momentum and site-decomposed DOS for Th 2 Al 
and Th 2 AlH 4 in Fig.|6]. DOS features for Th 2 Al and Th 2 AlH 4 show close similarity. Both 
exhibit metallic character since there is finite DOS at Ep. From the DOS histogram we 
see that Ep is systematically shifted towards higher energy in Th 2 AlH 4 . This is due to the 
increase in the number of valence electrons when Th 2 Al is hydrogenated. DOS for both 
Th 2 Al and Th 2 AlH 4 lie mainly in four energy regions (a) the lowest region around — 20 eV 
stems mainly from localized or tightly bound Th 6p states, (b) the region from —9.25 to 
—2.5 eV originates from bonding of H Is, Al 3p and Th 6d (Th-6<i and Al-3p states in Th 2 Al), 
(c) the region from —2.5 to OeV comes from bonding states of Al 3p and Th 6d and (d) the 
energy region just above Ep (0 to 3.5 eV) are dominated by unoccupied Th 4/ states. 

The semi-core Th 6p states are well localized and naturally their effect on bonding is 
very small. On comparing the Th 6p DOS of Th 2 Al and Th 2 AlH 4 , it is seen that the 
width is significantly reduced in Th 2 AlH 4 owing to the lattice expansion and the inclusion 
of additional energy levels below Ep. In the valence band region, the band width and 
DOS are larger for Th 2 AlH 4 than for Th 2 Al. Hydrogenation enhances interaction between 
neighboring atoms and thereby increases the overlap of orbitals and in turn results in the 
enlarged valence band width in Th 2 AlH 4 . In particular, the strong hybridization between 
Th 6d and H Is states increases the valence band width from 7.1 eV in Th 2 Al to 8.4 eV in 
Th 2 AlH 4 . Th 6d, Al 3p and H Is states are energetically degenerate in the valence band 
region indicating a possibility of covalent Th-H, Th-Al and Al-H bonds. However, the spatial 
separation Th-Al (3.22 A) and Al-H (3.02 A) is larger than Th-H (2.26 A). Therefore, covalent 
bonds between the former part is small whereas there is a significant covalent contribution 
between Th and H. In conformity with this the COHP and charge density analyses show 
directional bonding between Th and H (see Sec. |If IE2| and [1 f IE3|) . The accommodation 



of H in the interstitial position between Th and Al creates new bonding states between Th 
and H. This also enhances Th-Al distance around 2.2% compared with that in Th 2 Al. The 
consequence of this enhancement is that the Al DOS at the valence band region becomes 
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narrow and the splitting between the Al 3s and 3p states is almost doubled (see Fig.[|). The 
finite DOS at Ep which gives the metallic character of Th 2 Al and Th 2 AlH 4 comes from The/ 
states in addition to some states of Al p character. 

Another interesting feature of the total DOS of TI12AIH4 is presence of a deep valley 
around Ep which is termed as a pseudogap. Pseudogan features are formed not only in 
crystalline solidsll but occur also in amorphous phasesi! and quasicrystals.il Two mech- 
anisms have been proposed for occurrence of pseudo gap in binary alloys, one attributed 
to ionic features and the other to the effect of hybridization. Although the electronegativ- 
ity differences between Th, Al and H are noticeable, they are not large enough to explain 
the pseudogap in Th 2 AlH 4 . Hence hybridization must be the cause for creation of the 
pseudogap in Th 2 AlH 4 . There is a correlation between the occurrence of pseudogaps and 
structural stability^, that materials which possess pseudogaps in the vicinity of Ep usually 
have higher stability. This may be the reason for the higher value of AH in Th 2 AlH 4 than 
in Th 2 Al (Table fTJ). 



2. Charge density 

The analysis of the bonding between the constituents will give better understanding 
about the anisotropic changes in the structural parameters on hydrogenation of Th 2 Al. 
Fig. 0) shows the calculated valence charge density (obtained directly from the self-consistent 
calculation) within ab and ac planes for Th 2 AlH 4 . Th, Al and H atoms are confined to layers 
along c, Th and Al being situated in alternating metal layers with hydrogen in between, 
hence establishing a sequence of Th-H-Al-H-Th-H-Al-H-Th layers (see Fig. [I]). The H atoms 
are arranged in a chain like manner within the ab plane as also evident from Fig. 0b. It is 
interesting to note that the nature of H-H bonding is quite different along a and c. Although 
the H-H distance is 2.34Awithin the basal plane and 1.95 Aperpendicular to the basal plane, 
the bonding between the H atoms is not totally dominated by the latter interactions. In fact 
the CO HP analysis (sec. D 1 E 3| ) shows that the covalent H-H interaction within the ab plane 



is larger than within ac plane. The electron distribution between Al and H suggests ionic 
bonding between them, in line with their electronegativity difference of 0.7. In conformity 
with this the integrated charge inside the Al sphere is around 0.59 (0.8 according to the 
TBLMTO method) electrons larger in Th 2 AlH 4 than in Th 2 Al. 

The bonding between Th and H is predominantly covalent as evidenced by the finite 
charge between these atoms. The H-s electrons are tightly bound to the Th-d states, and 
Th-H arrangement forms a H-Th-H dumb-bell pattern. Now we will try to obtain a possible 
explanation for the short H-H distance within ac plane of Th 2 AlH 4 from the charge density 
analysis. The strong covalent interaction between Th and H in the ac plane(see Fig.|7|b) and 
the dumb-belled pattern tend to draw the electrons of H towards Th leaving only a small 
amount of electrons between the H along c to repul each other. The main reason for this 
short H-H distance is then reduced repulsion rather than bonding interaction between them. 



S 



3. COHP 



COHP is an extremely useful tool to analyze covalent bonding interaction between atoms, 
the simplest approach being to investigate complete COHP between the atoms concerned, 
taking all valence orbitals into account. COHP between Th-H, Th-Al, Al-H and H-H in 
TI12AIH4 are given in Fig.^. 

Owing to the very different interatomic distances between the H atoms in the ab and 
ac planes, special attention is paid to COHP in these planes. Both bonding and antibond- 
ing states are present almost equally in the VB region indicating that covalent interaction 
between the H atoms is not participating significantly to the stability of Th 2 AlH 4 . On the 
other hand, the bonding states are present in the whole VB region in COHP of Th-Al and 
Th-H indicating that covalent interaction between these pairs is contributing to structural 
stability. The presence of the large bonding states in the VB region of COHP for Th-H is 
the main reason for large value of heat of formation for Th 2 AlH 4 compared with Th 2 Al. In 
order to quantify the covalent interaction between constituents of Th 2 AlH 4 we have inte- 
grated the COHP curves up to E F for Th-Al, Th-H and Al-H giving -0.778, -1.244 and 
—0.072, respectively. Owing to the presence of both bonding and antibonding states below 
Ep in COHP the integrated value for H-H becomes negligibly small (—0.086 and —0.011 
within the ac and ab plane, respectively, but as the integrated value of bonding states alone 
is —0.571 and —0.136 respectively, the bonding H-H interaction is quite different in the 
two planes). Hence, one can conclude that the bond strength between the constituents of 
TI12AIH4 decrease in the order Th-H > Th-Al > Al-H > H-H. 

The experimental^ and theoretical studies show highly anisotropic changes in the lattice 
expansion on hydrogenation of Th 2 Al. According to the crystal structure of Th 2 Al the 
interatomic distance between the interstitial regions where one can accommodate H in the 
ab plane is 2.4 A. Hence, there is a large flexible space for accommodation of the H atoms 
in this plane without the need to expand the lattice. In contrast, the interatomic distance 
between the interstitial regions in the ac plane is only 1.65 A. So, large expansion of the 
lattice along c is necessary to accommodate H within the ac plane. As a result, even with 
a short H-H separation of 1.95 A, a volume expansion of 12.41% is needed when Th 2 AlH 4 is 
formed from Th 2 Al. The experimental observation of 0.105% lattice expansion along a and 
12.15% along c is found to be in excellent agreement with theoretically obtained value of 
0.03 and 12.41%, respectively. 



IV. CONCLUSION 

This study reports a detailed investigation on the electronic structure, bonding nature 
and ground state properties of Th 2 Al and Th 2 AlH 4 using first-principle method. The fol- 
lowing important conclusions are obtained. 

1) Th 2 Al and Th 2 AlH 4 are formed in the CuAl 2 -type crystal structure, the optimized atomic 
positions and lattice parameters are in very good agreement with recent experimental re- 
sults. 

2) Structural optimization gave a shortest H-H separation of 1.95 A, which is close to the 
recent experimental value of 1.97 A. 
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3) We observed a highly anisotropic volume expansion of 12.47% of the Th 2 Al matrix on 
hydrogenation to Th 2 AlH 4 , of which 99.76% volume expansion occurs perpendicular to the 
basal plane and neglible change along the basal plane. 

4) The large difference in interatomic distance between the interstitial regions within the 
ab and ac planes and the strong covalent interaction between Th and H along c keeps the 
H atoms close together in the c direction. This is the main reasons for highly anisotropic 
volume expansion on hydrogenation in TI12AI. 

5) Charge density and COHP analysis revealed that the Th-H bonds are stronger than the 
H-H bonds and other localized bonds in this structure. The formation of strongly bonded 
ThH 2 subunits in Th 2 AlH 4 makes repulsive interaction between the H atoms smaller along 
c and this is the precise reason for the violation of 2 Arule. 

6) There is a correlation between c/a and the structural stability of hydrated CuAl2-type 
phases. For phases with c/a < 0.825 the symmetry changes from I^/mcm to P 4 /ncc on hy- 
drogenation, whereas for c/a > 0.825 the crystal symmetry is not affected on hydrogenation. 

7) Density of states and bandstructure studies show that Th 2 Al and Th 2 AlH 4 are having 
non vanishing N(E^), resulting in metallic character. The cohesive energy analysis show 
that, Th 2 AlH 4 is more stable than Th 2 Al. 
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TABLES 



TABLE I. Atomic position of Th 2 Al and Th 2 AlH 4 





Th 2 Al 


Th 2 AlH 4 




X 


y 


z 


X 


y 


z 


Th Theory 
ExpH 
Expi 
Expi 


0.1583 
0.1588 


0.6583 
0.6588 


0.0000 
0.0000 


0.1632 

0.1656 
0.162 


0.6632 

0.6656 
0.662 


0.0000 

0.0000 
0.0000 


Al Theory 
Exp0 
Expi 
Expi 


0.0 


0.0 


0.25 


0.0 


0.0 


0.25 


0.0 


0.0 


0.25 


0.0 
0.0 


0.0 
0.0 


0.25 
0.25 


H Theory 
Expi 
Expi 




0.3705 
0.3707 
0.368 


0.8705 
0.8707 
0.868 


0.1512 
0.1512 
0.137 



TABLE II. Lattice parameters and inter atomic distances of Th 2 Al and Th 2 AlH4(all values 
are in A). 





Th 2 Al 




Th 2 AlH4 




Theory 


Exp, 


j3 


Theory 


Exp.H 


Exp.l 


a 


7.602 


7.618 


7.604 


7.626 


7.629 


c 


c = 5.723 


5.862 


6.433 


6.515 


6.517 


c/ a 


0.753 


0.769 


0.846 


0.854 


0.854 


Th-H 








2.273 


2.305 


2.387 


Th-Al 


3.199 


3.219 


3.269 


3.278 


3.291 


Th-Th 


3.403 


3.421 


3.509 


3.571 


3.495 


Al-H 








3.051 


3.061 


3.072 


Al-Al 


2.861 


2.931 


3.216 


3.257 


3.258 


H-H (ac-plane) 








1.945 


1.971 


1.790 


H-H (a6-plane) 








2.344 


2.305 


2.495 



TABLE III. Ground state properties of Th 2 Al and Th 2 AlH 4 



Compound 


Th 2 Al 


Th 2 AlH 4 


- AH ( in kJ mol" 1 ) 


18.35 


29.50 


E coh (eV/&tom) 


0.15 


0.185 


A r (^)(states/Ry-cell) 


58.42 


41.13 


So(GPa) 


93.42 


111.36 


B'o 


3.41 


3.48 
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FIGURES 



FIG. 1. The crystal structure of TI12AIH4. Five Th in face sharing tetrahedral configuration 
surrounding two hydrogen. Legends to the different kinds of atoms are given on the illustration. 

FIG. 2. Total energy (a) vs. c/a and (b) vs. unit cell volume for Th 2 Al where AE=E+ 106632. 

FIG. 3. The total energy (a) vs. c/a and (b) vs. unit cell volume for TI12AIH4 where 
AE=E+ 106636. 

FIG. 4. c/a for CuAFj-type phases and their corresponding hydrides. 
FIG. 5. Electronic band structure of (a) TI12AI and (b) TI12AIH4. The Fermi level is set to zero. 

FIG. 6. Total, site and orbital projected density of states for (a) TI12AI and (b) TI12AIH4. 

FIG. 7. Valence electron charge density plot for TI12AIH4 in the ab plane with 40 contours 
drawn between and 0.25 electrons/a.u. 3 

FIG. 8. Valence electron charge density plot between the Th and H atoms for TI12AH4 in the 
ac plane with 40 contours drawn between and 0.25 electrons/a.u 3 . 

FIG. 9. COHP of Th 2 AlH 4 , depicting the contributions from Th-Al, Th-H, Al-H and H-H 
interactions. The COHP for H atoms in the ab plane and ac plane are given in solid line and 
dotted lines respectively. 
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